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A direct generalization of the transverse momentum integrated(TMI) light-cone wave func- 
tion to define a transverse momentum dependent (TMD) light-cone wave function will cause 
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■ light-cone singularities and they spoil TMD factorization. We motivate a definition in which 

the light-cone singularities are regularized with non-light like Wilson lines. The defined TMD 
light-cone wave function has some interesting relations to the corresponding TMI one. When 
the transverse momentum is very large, the TMD light-cone wave function is determined per- 
turbatively in term of the TMI one. In the impact 6-space with a small 6, the TMD light-cone 
wave function can be factorized in terms of the TMI one. In this letter we study these relations. 
By-products of our study are the renormalization evolution of the TMI light-cone wave function 
\ and the Collins-Soper equation of the TMD light-cone wave function, the later will be useful for 

resumming Sudakov logarithms. 
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Exclusive B-decays play an important role for testing the standard model and seeking for new 
Q^, physics. Experimentally they are studied intensively. Theoretically, there are two approaches 
Q_ij of QCD factorization for studying these decays. One is based on the collinear factorization pQ, in 
which the transverse momenta of partons in a B-meson are integrated out and their effect at leading 
twist is neglected. The collinear factorization has been proposed for other exclusive processes for 
long time Another one is based on the transverse momentum dependent(TMD) factorization [2], 
^ ■ where one takes the transverse momenta of partons into account at leading twist by meaning of 
TMD light-cone wave function. The advantage of the TMD factorization is that it may eliminate 
end-point singularities in collinear factorization [I] and some higher- twist effects are included. The 
knowledge of the TMD light-cone wave function will provide a 3-dimensional picture of a S-meson 
bound state. However, it is not clear how to define the TMD light-cone wave function in a consistent 
way to perform a TMD factorization because of light-cone singularities 5 . 

In the collinear factorization the light-cone wave function for a 5-meson moving in the in- 
direction with the four velocity v is defined asjS]: 

<M* + ./*) = J ^-e lk+z ~ (0\q(z-n)LUoo,z-nh + ^L n (oo, 0)h(0)\B (v)) , (1) 

where we used the light-cone coordinate system, in which a vector is expressed as = 
(a + ,ar,a±) = ({a + a 3 )/V2,(a° - a 3 )/V2,a\a 2 ) and a\ = (a 1 ) 2 + (a 2 ) 2 . The field h{x) is 
the field for 6-quark in the heavy quark effective theory (HQET) and q{x) is the Dirac field for a 
light quark in QCD. The gauge link is L n defined with the light-cone vector = (0, 1,0,0) as: 

L n (oo, z) = Pexp ^— ig s J dXn ■ G(Xn + z)^j . (2) 
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In the definition the transverse momentum of q is integrated, resulting in that the field q and h are 
only separated in the light-cone direction. We will call $+(/c + ,^) as transverse momentum inte- 
grated(TMI) light-cone wave function. A direct generalization of Eq.(l) by undoing the integration 
to define the TMD light-cone wave function causes serious problems because it has light-cone sin- 
gularities, if quarks emit gluons carrying momenta which are vanishing small in the +- direction 
but large in other directions. Similar problems also appear in defining TMD parton distributions 
and fragmentation functions for inclusive processes. It has been shown that one can take in the 
definition the gauge link in the direction off the light-cone direction n and the factorization of 
inclusive processes can be done without light-cone sing ularities[3 El El- In this letter we propose 
to define the TMD light-cone wave function with gauge links slightly off the light-cone and study 
its relation to the TMI light-cone wave function. 

We introduce a vector = (u + ,u~ ,0,0) and define the TMD light-cone wave function in the 
limit u + « u~: 

Mk + ,k ± ,(,») = J ^^e^ + ^-^^M0|^)4(oo,z) 7 + 75^(oo,0)/ l (0)| J B( ? ;))|, + =0 , (3) 

the gauge link L u is defined by replacing the vector n with u in L n . This definition has not the 
mentioned light-cone singularity, but it has an extra dependence on the momentum k + through 
the variable £ 2 = 4(u ■ k) 2 /u 2 . This extra dependence is useful. The evolution in ( is controlled 
by the Collins-Soper equation 7_ which leads to the so-called CSS resummation formalismjB], and 
it will be derived here. The evolution with the renormalization scale [i is simple: 

<9c/> + + ,/c_L,C,/x) . . , + , . . 

V q- = (7ff+70)^+(* »*±»C»M)» ( 4 ) 

where 7 g and 7q is the anomalous dimension of the light quark field q and the heavy quark field h 
in the axial gauge u ■ G = 0, respectively. 

If one integrates kj_ in <j>j r {k + , k±, £, fx), the TMD light-cone wave function will in general not 
reduce to the TMI one. The reason is that the integration over k± in Eq.(l) has ultraviolet 
divergences, a renormalization is needed. Hence the integration over the transverse momentum 
in Eq.(l) is done in d — 2 dimension, if one uses (i-dimensional regularization, and then a UV 
subtraction is performed. In contrast, the transverse momentum k± in cj) + {k + , k±, £, /i) is in the 
physical space with d = 4, ultraviolet divergences will be generated if one integrates over k± and 
they are not subtracted in Eq.(3) because it is a distribution of k±. Therefore, one can not simply 
relate 4>+(k + , k±, £, jx) by integrating k± to the TMI light-cone wave function in Eq.(l). However, 
the TMD light-cone wave function has some interesting relations to the TMI one. If the transverse 
momentum k± carried by the parton q is much larger than the soft scale Aqcd, the 6-quark as a 
parton will also carry large transverse momentum because the momentum conservation. This can 
only happen if hard gluons are exchanged between the two partons and the exchange can be studied 
with perturbative QCD. Without the exchange of hard gluons, one expects that the partons will 
carry k± with a typical value at order of Aqcd- In the case with large k± the TMD light-cone 
wave function is determined in term of the TMI one as: 

Mk + ,k±,(,v)= dq + C ± (k + ,q + ,k ± ,()<5> + (q + ,v)+0((kir 2 ), (5) 
Jo 

where the function C± can be determined by perturbative QCD. By power counting^J C± is 
proportional to (A;^) -1 . When we consider the Fourier-transformed TMD light-cone wave function 
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into the impact 6-space: 

(p+(k + , b, £, (J>) can be related to the TMI one for small b as: 

Mk + ,b,(,n)= dq + C B (k + ,q + ,b,C,^ + (q + ^)+0(b), (7) 
Jo 

where Cb can also be calculated with perturbative QCD, i.e., it does not contain any soft divergence. 
Hence the relation represents a factorization. The leading order result is Cg = 5{k + — q + ). A 
similar factorization between TMD- and usual parton distribution was proven in '7,. We will 
determine the relation in Eq.(5) and Eq.(7) up to order of a s and show that the factorization in 
Eq.(7) holds at one-loop level. In determination of these relations we also derive the //-evolution 
equation of the TMI light-cone wave function and the Collins-Soper equation of the TMD one. 

The relation in Eq.(5) is useful for constructing models of the TMD light-cone wave function. 
The importance of the relation in the 6-space in Eq.(7) and the Collins-Soper equation is the 
following: When the TMD factorization is formulated in the 6-space, the relation allows to use the 
TMI light-cone wave function, while the Collins-Soper equation is used to resum large logarithms. 
Hence it is possible to have relations between two factorization approaches under certain conditions. 
It should be noted that other definitions of a TMD light-cone wave function are possible. A different 
definition is given in ^H] with a complicated structure of gauge links. With this definition the 
relations in Eq.(5) and Eq.(7) can also be studied. 

The functions C± and Cb are free from any soft divergence, i.e., infrared- and collinear singular- 
ities, we can use a partonic state instead of a B-meson state to determine them. We use the partonic 
state \b[mhV + kb),q(k q )), the momenta are given as k 1 ^ = {k q , k~, k q ±) and k£ = (k^ , k^, —k q ±). 
These partons are on-shell, i.e., k q = m q and v ■ fc& = in HQET. The variable k + of the 
wave functions is from to oo. Actually, from the momentum conservation, it is from to 
P+ = mi)V + + k^ + k q . Under the limit — > oo we have P + — > oo. If we set P + to be oo 
at the beginning, it will results in an ill-defined distribution like l/(fc + — fcj")+ for k + going to oo. 
Therefore we will take a finite P + in the calculation and take the limit P + — » oo in the final result. 
We illustrate this in detail for $ + . 





(d) 



Figure 1: Diagrams of one-loop contributions. Thick lines stand for 6-quark, double lines represent 
gauge links. 



At leading order & + (k + ,[i) = v(k g )"f + ^u v (ki ) )5(k + — fe+). At order of one loop, the contri- 
butions are from Feynman diagrams shown in Fig.l and Fig. 2., except from Fig. Id, Fig. 2b and 
Fig.2e. The contributions from Fig. 2. are proportional to the tree-level result. We use dimensional 
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regularization with d = 4 — e for U.V. divergence and give gluons a small mass A for infrared 
divergences. There are light-cone singularities in diagrams where the gluon is attached to a gauge 
link. They are cancelled between different diagrams. To show this, we take Fig.lc. as an example. 
The contribution from Fig.lc after the integration of the gluon momentum is: 



<M^)lic = - 2 ^v{k q ) 1 + 1 ,u{k b )e{k+ - fc+)(V) £/2 r(s/2) ^| ,!f r/2 > 



(8) 



where x = kg /k + and ^ 2 = 4(v k^) 2 . This contribution has the light-cone singularity at k + = kg . 
The contribution from Fig. 2c reads: 

$+(*:+ //)| 2c : 



2a, 



5{k+ - k+)v(k q ) 7 + l5 u(k b )^-(47r^ 2 r(s/2) 



dp' 



(x 2 + C 2 v (p + /K) 2 r £/2 



p 



(9) 

The integral is divergent for p + — > oo and p + — > 0. The divergence at p + = is the light-cone 
singularity and will be cancelled with that from Fig.lc. If we set P + to be oo at the beginning, the 
sum of these two contributions integrated with a test function f(k + ) is proportional to 



k+ 



kq 



In- 



kt) 2 



(10) 



This integral is divergent, because the divergence at p + — > oo in Eq.(9) is still there. The divergence 
has a geometrical reason jT^j. In HQET the 6-quark field h in Eq.(l) can be represented as a gauge 
link along the direction v and forms a Wilson line combining the gauge link L n . This Wilson line 
has a cusp singularity at the origin where the two gauge links ioin|14j. This divergence needs to 
be renormalized. With a finite P + the sum becomes proportional to the expression instead of the 
integral in Eq.(10): 



P+ 



k+-kt 



c 2 (k+ - k+f 



In 2 



C 2 (P" 



K) 2 



(ii) 



One can easily check that the expression is finite under the limit P + — > oo. Evaluating all diagrams 
we have the result at one loop with xp = P + /Ad~: 



$ + (k + ,n) 



v(k q )^ + ^f 5 u(k b )6(k 
2a 
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4 

K) x e v {i-x) 2 
fi 2 

the contribution from Fig. la is U.V. finite and is not needed for deriving our final results, as 
explained later. The above expression is a distribution for < k + < P + . From the above result 
one derives the renormalization evolution under the limit P + — > oo: 
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dq + 'j + (k + ,q + ,fj 1 )(j) + (q + ,Li), 



In- 



2v-k- 



-)5(k + -q+) + 



k + 9(q-< 



k+) + 9{k^ 



q+(q+ - k+) k 



(13) 
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In deriving this one should be careful with the plus description acting on different distribution 
variables. The plus prescription above is for q + . This result is in agreement with that in by 
noting the fact that the wave function defined in |12| is with the decay constant in HQET. From 
our explicit result it is observed that the wave function can not be normalized under the limit 
P+ —* oo as first observed in 6 and in other studies [TH]. 

To determine the mentioned relations we need to calculate the wave function <p + at one-loop 
order. All diagrams in Fig.l and Fig. 2 give contributions. The calculation is straightforward in the 
momentum space. The light-cone singularity is regularized by a small but finite u + . The detailed 
calculation and result will be presented elsewhere and we will only give final results here. To 
determine the function Cj_ we only need to calculate Fig. lb, Fig.lc and Fig. Id. The contribution 
from Fig. la is proportional to l/(fc^) 2 for large k± and will not contribute to C±. The function is 
determined by taking the limit of large k± and then P + — > oo. We obtain: 



' k+6{q + -k + ) 6{k + -q + ) 
q+{q+ - k+) + k+ -q+ 



+ 8(k^ 



k 



1) 



_L 



• (14) 









-4 


(a) 


(b) 


(c) 




-4 -j 




(d) 


(e) 


(f) 



Figure 2: The self-energy corrections. 



To determine the function Cb we do not need to calculate the contribution from Fig. la, again. 
The reason is that the contribution cj) + (k + , k±, (, fj)\ a is U.V. finite when integrated over k±. That 
means that cf) + (k + , b, (, fx)i a = 4>+{k + , b = 0, (, fj,)i a + 0(b) for 6^0. For small b we have: 



r 2ry r 

Mk + ,b,(,v) = v(k q h + 7 5 n(k b ){5(k + -k+) + -^[5(k + -k+)( 



3vr 

1 1 C 2 

-In 2 Cvb 2 (x P - l) 2 + 2 In fi 2 b 2 + - In \ 
4 2 [i 2 



m — 7T + - In — 7 



\n 2 C 2 b 2 --ln f , 2 b 2 ln 



m 



mi 



+ 



6(k^ 



A-1 



k+ -k+ 



ln(6 2 C(l - xY) - ln(6 2 m 2 (l - x) 2 ) 



q ^ ""q 

k + 9(k+ - k+) 



2--) 
3 ' 



k+{k+ -k+) J + 



)J} 



-<P + (k + ,b = 0,(,v)ia + O(b), 



(15) 
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with b 2 = 6 2 e 27 /4. This result is a distribution for < k + < P + . From the above result one can 
derive the Collins-Soper equation: 
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(-Q-:(l>+(k + ,b,(,Li) 



4a Sl CWT 1 2 «*i ^ e 

In in — 7T 

3vr 4 3vr ( 2 



<P+( k+ ,b,(,[j,). 



(16) 



The first factor is the famous factor K + (30 IE], the last factor comes because we used HQET for 
the heavy quark. 

Comparing the result in Eq.(12) and Eq.(15) and noting the fact that (p + (k + ,b = 0, CjZ-Ola is 
just n)ia, we can derive the function Cb under the limit P + — ► oo: 



C B (fc + ,g + ,&,C) 



5{k + - 
+8{k A 



3tt 
4 C 2 



£;+%+-£;+) , 9{k 



+ 



9 n 



1 



ln^(C 2 6 2 ) + ln(^6 7 ) + - In 



-fc+ 

M 2 



7vr2 
24 



;i7) 



which does not contain any soft divergence and does not depend on v. In the above the plus 
prescription is for q + . 

To summarize: We proposed the definition in Eq.(3) for the TMD light-cone wave function of a 
.B-meson. The definition does not contain the light-cone singularity and can be used for performing 
TMD factorization in a consistent way. Two relations between TMD- and TMI light-cone wave 
function are found. One is that the TMD light-cone wave function with large k± is determined 
by the TMI one. This relation will be important for constructing models of the TMD light-cone 
wave function. Another one is the factorization relation between the TMI light-cone wave function 
and TMD one in the impact b space with the small b. In studying these relations we also obtained 
the renormalization evolution of the TMI light-cone wave function and the Collins-Soper equation 
of the TMD one. The equation and the relation in the impact space are important. When TMD 
factorization is formulated in the impact space, the relation allows us to use the TMI light-cone 
wave function and the equation allows us to resum large logarithms. These issues will be discussed 
in another publication in the near future. 
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